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Abstract In 1665, Christiaan Huygens reported the observation of the synchroniza- 
tion of two pendulum clocks hanged on the wall of his workshop. After synchro- 
nization, the clocks swung exactly in the same frequency and 180" out of phase — 
anti-phase synchronization. Here, we propose and analyze a new interaction mech- 
anism between oscillators leading to exact anti-phase and in-phase synchronization 
of pendulum clocks, and we determine a sufficient condition for the existence of an 
exact anti-phase synchronization state. We show that exact anti-phase and in-phase 
synchronization states can coexist in phase space, and the periods of the synchro- 
nized states are different from the eigen-periods of the individual oscillators. We 
analyze the robustness of the system when the parameters of the individual pendu- 
lum clocks are varied, and we show numerically that exact anti-phase and in-phase 
synchronization states exist in systems of coupled oscillators with different param- 
eters. 



1 Introduction 

In 26 February 1665, Christiaan Huygens, in a letter to his father, fTj, reported the 
observation of the synchronization of two pendulum clocks closely hanged on the 
wall of his workshop. After synchronization, the clocks swung exactly in the same 
frequency and 180" out of phase. For attachment distances, less than 1 meter, the 
clocks always synchronize with the 180" phase difference. For larger attachment 
distances (> 4.5 m), synchronization, if it occurs, takes longer times. Huygens also 
noted that if the two clocks were hanged in such a way that the planes of oscillation 
of the two pendulums were perpendicular, then synchronization did not occur The 
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Huygens observations were the first time that synchronization effects have been 
described scientifically. 

Huygens justified the observed synchronization phenomena by the "sympathy 
that cannot be caused by anything other than the imperceptible stirring of the air 
due to the motion of the pendulum". 

In recent years, there has been a growing interest in the detailed analysis of syn- 
chronization phenomena, both from the theoretical and the experimental points of 
views. From the experimental point of view, Bennett et a/., |j2l, built an experimental 
device consisting of two interacting pendulum clocks hanged on a heavy support, 
and this support was mounted on a low-friction wheeled cart. This device moves 
by the action of the tensions due to the swing of the two pendulums, and the inter- 
action between the two clocks is caused by the mobility of the heavy base of the 
clocks. With this device, the anti-phase synchronization mode is reached when the 
difference between the natural or eigen-frequencies of the two clocks is less than 
0.0009 Hz. If the difference between these frequencies is larger than 0.0045 Hz, the 
two clocks do not synchronize, running "uncoupled" or in a state of beating death, 
||2|. This situation is unsatisfactory when compared with the observations of Huy- 
gens. For example, a difference of order of Ad) — 0.0009 Hz for the two pendulum 
eigen-frequencies corresponds to a difference in the lengths of the pendulum rods 
of the order of At = ^l^l^Aco/A-n, which gives, for ^ = 1 m and g = 9.8 ms^^, 
A£ = 4 mm, and for £ = 0.178 m (the length of the pendulum rods used by Huy- 
gens, (|T|), A£ = 0.02 mm, a precision that Huygens certainly could not achieve. 
According to Bennett et al. [2 p. 578], Huygens's results depended on both talent 
and luck. 

In the experiment of Bennett et al., 0, the in-phase synchronization is the natu- 
ral way of synchronization of the two pendulum clocks. However, due to a detailed 
description of the Huygens findings, we can believe that the in-phase synchroniza- 
tion was never observed by Huygens. 

Another experimental model mimicking the Huygens's clocks system, consists 
of two pendulums whose suspension rods are connected by a weak string, and one 
of the two pendulums is driven by an external rotor, |3| and |4|. In this system, 
the in-phase synchronization is approximately achieved with a small phase shift, 
and the experimental measurements and the model analysis both agree. The nu- 
merical results of Fradkov and Andrievsky for this device, |4|, show simultaneous 
and approximate in-phase and anti-phase synchronization, tuned by different initial 
conditions. In another experimental device made of two rotors controlled by ex- 
ternal torques (||5]|6l), Andrievsky et ah, (|5 |, reported approximate anti-phase and 
in-phase synchronization of the two oscillators. In this experiment, the synchroniza- 
tion parameter is the stiffness of a string connecting the two rotors. 

For demonstration purposes, Pantaleone, |7|, reported an experimental device 
constructed with two metronomes on a freely moving base. In this metronomes ex- 
periment, the phase difference of the synchronous state is close to 0". Increasing 
the damping effect on the freely moving base, the author also reported approximate 
synchronization with a difference in phase close to 180". From the theoretical point 
of view, the equations describing this experimental device lead to the Kuramoto syn- 
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chronization model, [8 |, where the synchronization mechanism is associated with a 
non-linear effect associated with the phase difference between the oscillators. 

In the Bennett et ai, fT\, and the Pantaleone, |7|, experimental systems, the in- 
teraction mechanism between oscillators is obtained by a moving base, an idea ad- 
vanced by Kortweg in 1906, |9|. 

In the experimental systems described above, there is no clear evidence of what 
mechanisms are in the origin of the anti-phase synchronization, as described by 
Huygens. In some experiments, it appears that if the pendulums have slightly differ- 
ent periods, the two oscillators may not synchronize, ||2JB1- 

However, as this special type of collective rhythmicity occurs in biological sys- 
tems and several other natural phenomena, [101, where individual periods are differ- 
ent, it is important to derive and to understand the interaction mechanisms leading 
to exact synchrony. Besides all the attention of the scientific community for this 
synchronization phenomenon, there is no clear evidence of a mechanism, leading to 
anti-phase synchronization. 

In this paper, we propose and analyze a new interaction mechanism between os- 
cillators leading to exact anti-phase and in-phase synchronization. The synchroniza- 
tion mechanism is obtained with a damped elastic string, and the oscillators under 
analysis can be simple harmonic oscillators, pendulums, and any type of non-linear 
oscillators with a limit cycle in phase space, as is the case of pendulum clocks, ifTTI . 
The main result of this paper is to show that exact anti-phase synchronization can 
always be achieved for systems of coupled oscillators. 

This paper is organized as follows. In section [2] we introduce the synchroniza- 
tion model in its full generality, discussing its physical assumptions, and we derive 
the equations of motion of the interacting oscillators. Then, we make the approxi- 
mation of small oscillations. In section|3] we introduce a simplified pendulum clock 
model, and we prove that the ordinary differential equation describing the dynamics 
of the clock has a limit cycle in phase space. Based on the linear model derived in 
section |2] in section]?] we discuss the concept of anti -phase synchronization, and 
we find a sufficient condition for the existence of an exact anti-phase synchronized 
state for the two pendulum clock system. This result is stated as Theorem[T] the main 
result of this paper After exploring numerically the phase space structure of the so- 
lutions of the model equations as a function of a control parameter, we show the 
existence of exact anti-phase and in-phase synchronization states for the Huygens's 
two pendulum clocks system. The tuning between the two types of synchronization 
regimes can be controlled through a damping constant associated with the (elastic) 
interaction mechanism and by the choice of initial conditions. These results justify 
some numerical results published previously, IITZl . In section |5] we analyze numer- 
ically the persistence of in-phase and anti-phase synchronization phenomena when 
we vary the parameters of the individual pendulum clocks. We show that the anti- 
phase and the in-phase synchronization regimes persist even if the two oscillators 
have different eigen-periods and parameters. Finally, in section |6] we resume the 
conclusions of the paper. 



4 



Rui Dilao 



2 A model for the synchronization of the two pendulum clocks 

In the Huygens two pendulum clocks system, the pendulums are hanged in a com- 
mon support, and the only possible interaction is due to the tension forces generated 
by the oscillatory motion of the two pendulums. These tension forces propagate 
through the common support, that we consider to be elastic. The role of the tension 
forces in the synchronization mechanism is corroborated by the Huygens's finding 
that when the planes of oscillation of the two hanged pendulums are mutually per- 
pendicular, no synchronization is observed. In fact, the components of the tension 
forces generated by the motion of the pendulums are in the plane of motion of the 
pendulums, and force the motion of the two attachment points. To be more specific, 
we consider the geometric arrangement of Fig. [T] where the two pendulums have 
masses mi and 1112, and lengths £1 and £2, respectively. 



Fig. 1 Model to analyze the synchronization of the Huygens's two-pendulum clocks system. The 
two pendulums are a representation of two nonlinear oscillators. The interaction between the pen- 
dulums is done by the tension forces at the attachment points, and they actuate through an elastic 
and resistive media. Each attachment points is considered to have mass M. 

The pendulums are considered connected by a massless string with stiffness con- 
stant k. The perturbations that propagate along the string are damped, and the damp- 
ing force is proportional to the velocity of the attachment points of the string, with 
damping constant p. The damped string simulates the interaction effects that propa- 
gate through the elastic and resistive media. We consider that the attachment points 
of the pendulums, located at the horizontal coordinates x~xi and x — X2, have equal 
masses, and we denote this mass by M. As we shall see, the introduction of this mass 
is necessary to obtain explicitly the equations of motion. 

The system of Fig.[T] considered without the damping forces, is described by the 
four degrees of freedom Lagrangian, 




mi 



m2 



L = \m\ {£\el +x\ + 2£ixi 0i cos 0i ) + mig£\ cos Oi 

+ \m2 {£l 0| +xj + 2^2^202 cos 02) + m2g£2 COS 02 

+ ^M{x\+xl)- ^k{x2-xi)^ , 



(1) 
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where 9i and 62 are the angular coordinates of the two pendulums, g is the acceler- 
ation due to the gravity force, and the last two terms describe the interaction mech- 
anism between the two pendulums. From ([T]), the Lagrange equations of motion of 
the two interacting pendulums are. 



mili Oi +migsin0i = — wii'i cos0i 
2 02 + m2g sin 02 = -^2x2 cos 02 
[M + m\)xi + mi ^101 cos Oi = m\l\ 0j^sin0i +k{x2 —x{) 

[M + m2)i2 + 'W2^2 02 cos 02 = '"2^2 02 ^2 ^ k{x2 — X\) . 

Introducing the damping effects into system of equations (|2]), we obtain, 

mii\d\ + /i(0i, 0i) +mi^sin0i — — mii'i cos d\ 

m2^2 02 + /2 ( 02 , 02 ) + mg sin 02 = -m2X2 cos 02 

{M + m{)x\ +2px\ +mi^i0icos0i = mii?i 0^^ sin 0i -\-k{x2—xi) 

(M + m2)x2 + 2px2 + 1112^2 02 cos 02 = m2^2 02 sin 02 — A;(x2 — JCi ) , 



(2) 



(3) 



where p is the damping constant of the attachment points of the pendulums, and the 
functions /i (0i , 0i ) and /2(02, 02) describe the escaping mechanism of the pendu- 
lum clocks (§|3]l. 

The system of equations Q implicitly defines a system of ordinary differential 
equations. If M > 0, > 0, m\ > 0, ^2 > and 1112 > 0, the system (|3]l can be solved 
algebraically in order to the higher derivatives, and we obtain. 



mit\d\ + /i(0i, 0i) +migsin0i = —m\ cos 0\F\ 

m2(.2 02 + /2 ( 02 , 02 ) + m2g sin 02 = -m2 cos 02^2 

xi = Fi 
X2 = F2 , 



(4) 



where, 

/i (01, 01 ) cos 01 +migsin0i cos 0i —2px\ +mi^i0[^sin0i +^(x2 — xi) 

^1 = ^ 9 

M + misin 01 

P _ /2(02,02)cos02 + m2gsin02COS02-2pi2 + m2402 sin02-fe(x2-xi) 

M + m2 sin^ 02 

(5) 

The system of ordinary differential equations (jSll are a synchronization model for 
the Huygens's two pendulum clocks system, lfT2f 7 Here, we will consider only the 
case M > 0. The case M = will be analyzed elsewhere. 

For small amplitude of oscillations in 0i and 02, the system of equations Q-Q 
simplify, and we obtain. 
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mi£l0l+/l(0i,0i)+mig0i = -^(/i(0i,0l)+mig0i-2pii+fc(x2-xi)) 

JU2. 

OT2^202+/2(02,02)+'«2g02 = {.fliOl, Ol) + migOl - 2pX2 - k{x2 - Xi)) 

JCi = -2pii+fc(x2-xi)) 

JC2 = ^ (/2(02, 62) +'W2g02 -2pi2 -fc(x2 -Xi)) . 

(6) 

In this paper, our goal is to analyze the synchronization properties of the solutions 
of the system of ordinary differential equations (j6]). 

To model the Huygens's two pendulum clocks experiment, we have to choose a 
specific form for the functions 0i) and /2(02, 62), describing the oscillatory 

behavior of pendulum clocks. 



3 A simple clock model 

To restore the energy lost by a pendulum clock during one period, the sustained os- 
cillations can be maintained by an impulsive force acting on the pendulum rod, or 
by the force originated by the smooth unwinding of a circular string attached to the 
pendulum balance wheel, ifTTl pp. 169-200]. In any case, the dynamics of a pendu- 
lum clock is modeled by a two-dimensional dynamical system with a limit cycle in 
phase space. This same qualitative behavior can be obtained with an oscillator actu- 
ated by a non-linear damping force, piecewise proportional to the angular velocity 
of the pendulum. For small amplitude of oscillations, the proportionality constant 
is positive, and, for large amplitudes of oscillations, the proportionality constant is 
negative. 

To simplify our analysis, we take, as a qualitative model for a pendulum clock, 
the following second order differential equation, 

mW+f{e;X,e)e+mge = o, (?) 

where. 

The function —f{9;X, 0)0 is the damping force of the pendulum clock, and X and 
are positive constants. 

Proposition 1. IfX > 0, > 0, m > 0, ^ > 0, and g>Q, the second order differential 
equation with the damping function ([S]), has a unique and stable limit cycle in 
phase space. 



Proof. Assume that A > 0, > 0, m > 0, ^ > 0, and g>Q. Define the function. 
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where, T] = A /{mi), and ap- — g/t. With the new coordinate, x~ 9 + F{9;ri ,9), 
the differential equation (|7]l, with damping function ([8]), has the Lienard form, 

= x-F{9;ri,9) 
= -(o'-9. 

where F{9;ri,0) is a continuous and odd function of 9. The existence and stabiHty 
of a limit cycle solution for equation (|9]l follows from the Lienard theorem, |l3l pp. 
179-181]. □ 

In the conditions of Proposition [T| the differential equation (|7]i, with damping 
function (|8]l, has a unique limit cycle in phase space. In Fig. |2] we show the limit 
cycle in phase space of this non-linear oscillator, for two values of the parameter 



a) b) 




-1.0 -0.5 0.0 0.5 1.0 -0.4 -0.2 0.0 0.2 0.4 

9 9 

Fig. 2 Limit cycles solutions of the differential equation with the damping function jsj. In a) 
we have chosen (O^ = 9.8, and in b), co^ = 0.005, where CO- = g/£, m = 1 and £= I. The parameter 
values of the damping function are, A = 0.1 and 6 = 0.1. The two vertical lines for 6 = ±0 show 
the discontinuity of the tangential derivative along the limit cycles. 



The non-linear oscillator model defined by the differential equation (j7]i, with 
damping function (|8]l, has all the qualitative properties found in more accurate pen- 
dulum clock models. 111], and will be used in the next sections to model Huygens's 
clocks. 
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4 Synchronization of two pendulum clocks with equal 
parameters 



By (j6]l, the two interacting pendulum clocks with equal parameters are modeled the 
system of differential equations (|6]), 



Ml 



(X2-Xl) 



1 111 o Jc 



(10) 



1 ~ ■ fn O 

,2- -f{e2;X,e)e2- -892+2^x2 



M 



{X2-Xl), 



where f{9;X,9) is given by and co^- = g/£. 

Our goals here is to show that the solutions of the system of differential equations 
( [T0| synchronize in anti-phase. For that, we begin by analyzing the solutions of the 
system of equations ( 10 1 for the particular case where the initial conditions in 9i and 
02 are bounded by 9, that is, 1 0i (0) | <9 and 1 62(0) | < 6. As the system of equations 
( lOi, together with Jsll, is piecewise linear, and if |0i(f)| < 9, and \92{t)\ < 9, for 



every t e [0,f*], where t* is a positive constant, then the piecewise linear system of 
equations ( [T0| can be written as the linear first order system of differential equations. 
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(11) 



where we have introduced the new variables, vi = 9\, V2 = 92, w\ = ii, 1^2 = i2, 
and the new constants, A = -0)^(1 + m/M), and B = 2A(1/(M^) + l/(m£)). As, 
|0i(f)| < and 1 02(f) I < 0' for every t e [0,f*], then, for the same initial conditions. 



the small amplitude solutions of ( |T(j] i and (111 coincide in the interval [0, t*]. 

Denoting by Q the matrix in we have, DetQ = 0. A simple inspection of 
Q, shows that Q has only one zero eigenvalue, corresponding to the eigendirection 
defined by the equation x\ = We investigate now the stability of the line of fixed 
points x\ = X2 of both equations ([TO]) and (111. 



Proposition 1. If X > 0, M > Q, t > Q, m > Q, and p > is sufficiently small, then 
the systems differential equations ( jiOp and have a line of fixed points with 



The vector = (0, 0, 0, 0, 1 , 0, 1 , 0) is such that Mx = 0. 
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coordinates, 9\ — 0, v'l =0, 02 = 0, V2 — Q,w\— 0, wi — 0,x\— X2, and this line 
of fixed points is Lyapunov unstable. 

Proof. We assume that, A >0, M>0, ^>0 and m > 0. The existence of the Hne of 
fixed points follows by solving the equation Qy — 0, where, 



y 



(01,Vi,02,V2,Xi,Wi,X2,W2). 



For p = 0, the characteristic polynomial of the matrix 2 in < 1 1 ' is, 

qp^Q{x)=x^{A + {B-x)x) {2gkm + M{A + {B-x)x) {Mx^ + 2k)t-AkxX) , 

where A = -0}^{\+m/M) and B = 2X{\/{Mt) + l/{m()). As the polynomial 
qp=o{^) has the two roots. 



with positive real parts, B > 0, then, by continuity of Q{p), for sufficiently small 
p > 0, the characteristic polynomial qp>Q{x) of the matrix Q, has also eigenvalues 
with positive real parts. Therefore, for sufficiently small p > 0, A > 0, M > 0, ^ > 0, 
m > 0, the line of fixed points of the linear system ( 11 1 is Lyapunov unstable. As 
both systems of equations 



lOi and (11 



have the same phase space orbits near the 
common line of fixed points, the local properties of the flows are the same, and 
this is sufficient to prove the local instability of the flow defined by the system of 
equations ( [TO) i. □ 

Proposition [2] gives the conditions of nonconvergence to the quiescent state of 
the solutions of the system of linear equation ([TTJ. This quiescent state is the line of 
fixed points xi X2 in the eight-dimensional phase space. Under the conditions of 
the Proposition|2] the two pendulum clocks have sustained oscillations, in the sense 
that, lim,^«, 0i (f ) ^ 0, and Hm,^oo 02(0 ^ 0. 

To find a sufficient condition of existence of anti-phase synchronization of the 
two non-linear oscillators, we assume now that the asymptotic solutions of the sys- 
tem of equations ( fTO] ) synchronize in anti-phase. As the two pendulum clocks are 
identical, with, ?„ = fo + nh, where n = 0, 1 , • • •, and fo is the initial time, we must 
have. 



hm 01 (f„) = - lim 02 (f„) 



Um x\{t„) 



lim X2{t„) 



hm vi(f„) = - hm V2(f„) 
hm w\{t„) = — hm W2(f„). 



(12) 



for every /i > 



^ In this context, we say that the two pendulum clocks synchronize in anti-phase if, for every h > 0, 
lim,„^„ 9[ {tn) = — lim,j,^co 02 (fn)- As we shall see in § |5] this definition can only be used in the 



case of identical pendulum clocks. 
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We now define the new variables, 9 — 9i + 62, v = vi + V2, x = xi + X2 and 
w = wi + W2- Then, adding the corresponding variables in the system of equations 
( fTT) , we obtain. 



V 
X 



( 
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M 
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2A 












2p 
M£ 
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_2p 



M j 



V 
X 

\wj 



(13) 



In the following, we call the system of equations ( 13 1 the reduced system of equa- 



tions associated to system (111. The reduced system of equations ( 13 1 has a line of 
fixed points with coordinates, 0==O, v = 0, w = and x ^constant. If the line of 
fixed points of the reduced linear system ( [T3] l is asymptotically stable, then we have. 



lim 0(f) = lim + lim 02(0 = 
limv(f) = limvi(f) + limv2(f) = 

t^oa r^oo /^oo 

lim w{t) — lim M'i(f) + lim W2{t) — . 



(14) 



If, for every h > 0, the first condition in ( 12 1 is verified, the first condition in ( 14 1 



is also verified, and synchronous solutions of equation ( 10 1 are also stable solutions 



of equations ( 13 1. 

As, by Proposition |2] lim,,j^oo 0i(f„) ^ 0, and lim,„^oo 02(fn) 7^ 0, for any initial 
condition away from the line of fixed points, any solution of the differential equation 



(10 1 that anti-phase synchronizes is also an asymptotically stable solution of the 
reduced system of equations 



13 1 



13 1 is asymptotically 
is unstable 
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Hence, if the line of fixed points of the reduced system 
stable and the line of fixed points of the linear system of equations 
(Proposition |2j, then the solutions 0i(f) and 02(f) of the system of equations ( 10 1 
anti-phase synchronize. The asymptotic stability condition for the reduced the linear 
system ( [T3] l, together with the instability condition of Proposition |2] both give a 
sufficient condition for the existence of exact anti-phase synchronization of the two 
identical non-linear pendulum clocks. 

To analyze the stability properties of the line of fixed points of the system of 
equations ( [T3] l, we calculate the characteristic polynomial of the matrix P in ( [T3] l, 



P{y) = y{mMly^ 
+y{m 

= yp\(y) 



f {2mep - 2mX - 2MX ) 
niMiop- - 4A p ) + 2mgp ) 



(15) 



To the eigenvalue x = corresponds the eigendirection e^. This eigendirection is the 



line of fixed points of the linear system ( 13 1. If all the eigenvalues of the polynomial 



p\{y) in ( 15 I have negative real parts, any initial condition away from the line of 
fixed points x =constant, evolve in time to this line of fixed points. Therefore, we 
have: 
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Theorem 1. We consider the system of differential equations [lO), with damping 



function If X > 0, 9 > 0, m > 0, £ > 0, M > 0, k > 0, g > 0, and p > is 



sufficiently small, and if the reduced linear differential equation ( 13 1 has only non- 
positive eigenvalues, then, the solutions of equation ( jiOp , with damping function ([S]), 
synchronize in anti-phase, in the sense that, for every h > 0, 

lim 01 (f„) = - lim 02 (f«), 

where, tn = fo +nh, lim,,j^oo 0i (?«) 7^ 0, and lim,jj^oo 02(fn) 7^ 0. Moreover, if the 
polynomial, 

q{p) = 4-mlXp^-{4-mX^ + 4-MX^ + gm^i)p 
+gm^X + Igm^MX + gmM^X , 

has two real roots pi and P2, and p obeys to the inequalities, 

P\<p <P2 

P>Po=f(l + ^), 

then, for any initial condition away from the line of fixed points of the system of 
equations {10), the solutions of equation {10\, with damping function ([S]), anti-phase 
synchronize. 

Proof. The sufficient condition for the existence of anti-phase synchronization of 
the two pendulum clocks has been derived before the statement of the theorem. The 



instability of the line of fixed points of the system of equations (10 1 has been proven 
in Proposition |2] To prove the condition of non-positivity of the eigenvalues of the 
matrix in the system of equations ([T3]l, we use the Routh-Hurwitz criterion, 1 14|. By 



(1 5 1, as /? 1 (y) = a^y^ + a + a2y + 03 , and as, by hypothesis, oq > and 03 > 0, by 
the Routh-Hurwitz criterion, if ai > and {a\a2 — floaa) > 0, then the polynomial 
p\{y) has only roots with negative real parts. As, 

(fllfl2-fl0fl3) = ~q{p) = ~'XmtXp^ + {AmX'^ + AMX'^ + gm^l)p 
—gm^X — Igm^MX — gmM^X , 

the polynomial q{p) has a global minimum for positive values of p and can have 
two real roots. This proves the first inequality of the theorem. The second inequality 
follows from the Routh-Hurwitz condition, a\ — 2m£p — 2mX — 2MX > 0. □ 

Theorem [T] gives a sufficient condition for the existence of exact anti-phase syn- 
chronization in the Huygens's two pendulum clocks system. To test numerically the 
results of Theorem [T] we take for the parameters of the nonlinear oscillator (|7]l-(j8]l 
the values, g = 9.8, m = 1, ^ = 1, A = 0.1, and = 0.1. For these parameter values, 
the period of the solutions on the limit cycle of the equation (j7]i is, T = 2.008. The 
parameters describing the interaction between the pendulum clocks have been set to 
the values, k= 10, M 0.1, and p is a free parameter 

To test the conditions of Theorem [T] in Fig. |3] we have plotted the eigenvalues 



with the largest real part of the matrices Q and P, of the linear systems (111 and ( 13 1 
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respectively, as a function of the damping parameter p . The zero eigenvalue has been 
excluded from the characteristic polynomials of the matrices Q and P. According 
to Theorem[T] if, 0.121 = pi < p < p2 = 24.489, p > po ^ 0.11, and the matrix 
Q has eigenvalues with positive real parts, the two pendulum clocks synchronize 
in anti-phase. Numerically, for our reference parameter values, if p < p3 = 0.393, 
the matrix Q of the system of equations ( [TT] i has positive eigenvalues. Therefore, 
the conditions of Theorem [T| imply that, if p e [pi,p3], the two pendulum clocks 
synchronize in anti -phase, Fig.|3] 
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Fig. 3 Eigenvalues with the largest real part of the matrices Q and P of the linear systems and 
^13^ , respectively, as a function of the damping parameter p . The other parameter have been fixed 
to the values: g = 9.8, m = 1, ^ = 1, A = 0.1, 6 = 0.1, it = 10 and M = 0.1. If p > p, =0.121 and 
p < P2 = 24.489, all the eigenvalues of the matrix P have non-positive real parts, and the line of 
fixed points of system \13) is Lyapunov stable. If p < p3 = 0.393, the matrix Q has eigenvalues with 
positive real parts, and the line of fixed points of system |T0| is Lyapunov unstable. By Theorem[T] 
if p e [pi ,p3], the two pendulum clocks synchronize in anti-phase. 



By Theorem [T] the anti -phase synchronization between the two non-linear pen- 
dulum oscillators exists for p e [pi = 0.121,p3 = 0.393]. In Fig. |4j we show the 
time evolution of the angular coordinates and attachment points of the two pen- 
dulum clocks, for p = 0.37, and calculated numerically from the system of equa- 
tion ( 10 1. We have chosen for initial conditions the coordinate values, 0i(O) — 0.2, 
02(0) =0.3, XI (0) =0, X2(0) =0, 01 (0) =0, 02(0) = 0, ii(0) =0 andi2(0) = 0. 
The two non-linear pendulums and the attachment points synchronize in anti-phase. 
Numerically, the period of the synchronized state is T = 2.457, which differs from 
the eigen-period of the uncoupled non-linear oscillators, T = 2.008. The anti-phase 



synchronized state of the system of equations ( 10 1 corresponds to a stable limit cycle 
in the eight-dimensional phase space. 
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Fig. 4 Numerical solutions of the system of equations jlo| , with damping function ([8j, describing 
the coupling of two identical pendulum clocks. The parameter values of the simulation are: m = 1, 
£ = 1, g = 9.8, A: = 10, M = 0.1, X = 0.1, 6 = 0.1 and p = 0.37. The initial conditions are: 
01 (0) = 0.2, 62(0) = 0.3, xi(0) = 0, X2(0) = 0, 61(0) = 0, 62(0) = 0, ii(0) = andi2(0) = 0. 
In a) and b), we show the time evolution of the angular coordinates of the two pendulum clocks, 
before and after anti-phase synchronization, respectively. In c) and d), we show the asymptotic 
solutions in the reduced phase space of the two pendulum clocks (c), and of the two attachment 
points (d). For comparison, in c), we show the limit cycle solution (thin line curve) of the reference 
equation j7]l, with damping function jSj.The period of the anti-phase oscillations is T = 2.457. 
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Fig. 5 Anti-phase synchronization of the two pendulum clock system, for the same parameter 
values of Fig. [4] except for the damping parameter p that, in this case, has the value p = 0.1. In 
this simulation, the period of the exact anti-phase oscillations is T = 2.463. In c), the thin line 
curve is the limit cycle solution of the reference equation |7]l, with damping function |8]l. 
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Fig. 6 In-phase synchronization of the two pendulum clock system, for the same parameter values 
of Fig. [4] except for the damping parameter p that, in this case, has the value p = 0.01. In this 
simulation, the initial conditions are : 9, (0) = 0.2, 02(0) = 0.3, xi (0) = 0, ^2(0) = 0, 9, (0) = 0, 
02 (0) = 0, xi (0) = and X2 (0) = 0. The period of the in-phase oscillations is 7 = 0.606. In c), the 
thin line curve is the limit cycle solution of the reference equation with damping function ([8}. 
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Fig. 7 Anti-phase synchronization of the two pendulum clock system, for the same parameter 
values of Fig.|6] and p = 0.01. In this simulation, the initial conditions are : 61 (0) = 0.2, 62(0) = 
-0.3, xi{0) = 0, X2(0) = 0, ei(0) = 0, 62(0) = 0, ii(0) = and X2(0) = 0. Initially, the two 
pendulums are approximately in anti-phase. The period of the anti-phase oscillations isT = 2.463. 
In c), the thin line curve is the limit cycle solution of the reference equation (|7j, with damping 
function Isl. 
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Fig. 8 Coexistence of anti-phase and in-phase synchronization and aperiodic regimes for the 
asymptotic solutions of the system of equations {10\ . We show the limit cycle solutions of the sys- 
tem equations ^10^ for the same parameter values of Fig. [4] and several values of p: a) p = 0.075, 
b) p = 0.07, c) p = 0.06, and d) p = 0.02. We have plotted the asymptotic solutions in phase 
space for two different initial conditions. In one case we have taken 61 (0) = 0.2 and 61 (0) = 0.3, 
with the other initial conditions equal to zero. In the other case, we have taken 0i (0) = 0.2 and 
6\ (0) = —0.3. For p > 0.075, different initial conditions leads always to anti-phase synchroniza- 
tion. For p < 0.02, we have two limit cycles in phase space, one corresponding to anti-phase 
synchronization, and the other to in-phase synchronization. The thin line curve is the limit cycle 
solution of the reference equation jTJ, with damping function jSj. 



In Fig.|5] we have decreased the parameter p to values below the transition values 
Po and pi of Theorem[r| In this case, the two non-linear oscillators also anti-phase 
synchronize. The numerical integration for several initial conditions shows that, for 
the damping parameter value p =0.1, the system of equation ( 10 1 has a stable limit 
cycle in the eight-dimensional phase space. In this case, we are outside the condi- 
tions of Theorem[T] 

Decreasing further the parameter p, for p < 0.06, and for the same initial con- 
ditions as in Fig|4] the two oscillators synchronize with the same phase (in-phase). 
Fig. [6] However, changing the initial conditions for an approximate anti-phase ini- 
tial state of the two pendulums, we obtain anti-phase synchronization, Fig.|7] This 
simple fact shows that, for p < 0.07, there are two stable limit cycles in the eight- 
dimensional phase space, and these limit cycles have their own basins of attraction. 
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To analyze the transition from the anti-phase to the in-phase synchronization 
asymptotic regimes, we have chosen p in the interval [0.06,0.07], and we have 
changed the initial conditions of the numerical simulations. For p in the interval 
[0.06,0.07], the initial conditions that eventually lead to an in-phase synchronized 
regime show quasi-periodic behavior in time. If p < 0.06, the asymptotic in-phase 
synchronized state is a limit cycle in phase space. This suggests the existence of a 
non-local bifurcation for p near the value p ~ 0.07. 

In Fig. [8] we show the transition from the anti-phase to the in-phase synchro- 
nized asymptotic state as well as the bifurcation behavior of the attractors. We have 
superimposed in the same figure the two stable limit cycles that can be reached from 
different initial conditions. In Fig.lHh), p = 0.075, there is only one asymptotic anti- 
phase synchronized state. In Fig.pp), p = 0.020, there are two stable asymptotic 
synchronized state, one in-phase and the other in anti-phase. In Fig.jS]?) and[8};), we 
have a stable limit cycle for the asymptotic anti-phase synchronized state, but the 
asymptotic in-phase synchronized sate shows quasi-period behavior 

Therefore, we have shown that the interaction mechanism proposed in section |2] 
leads asymptotically in time to exact anti-phase synchronization of coupled identi- 
cal oscillators. For the same parameter values, we can have anti-phase and in-phase 
synchronization, determined by different initial conditions, and the tuning between 
the two types of synchronization regimes can also be controlled through the damp- 
ing constant p . Further numerical tests show that these results are still true in the 
singular limit M ^ 0. 



5 Synchronization of two pendulum clocks with different 
parameters: Robustness 

In the previous analysis, we have consider that both oscillators are characterized by 
the same parameters. However, in real experiments this is not realistic and we must 
consider the persistence of synchronization when the parameters of the pendulums 
are different. Here, the persistence of the anti-phase and the in-phase synchroniza- 
tion states is analyzed numerically. We take, mi — m, m2 — in{l + e), £i — £ and 
£2 — ^(1 + 5), where £ and 8 can have positive or negative values. In this case, the 
equations of motion (|6]l are rewritten as. 




k 

m 




(16) 



{x2-xi), 
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where A, 0) is given by jsjl. 
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Fig. 9 Anti-phase synchronization of two pendulum clocks with different lengths and masses. In 
a), we have the anti-phase synchronized state as in Fig. [5] In b), we show the anti-phase synchro- 
nized state for £ = 0.4 and 5 = 0.4. In c) and d), we show the limit cycles in phase space (thick 
lines) of the angular coordinates and of the attachment points of the two pendulums. The thin lines 
are the limit cycles for the cases e = and 5 = 0. In b), the period of oscillation isT = 2.661, and 
in a), T = 2.463. When the two non-linear pendulums have different parameters, the anti-phase 
synchronization condition, \im,^^a:, di{t„) = — lim(^^oo 02(?n), derived in §|4] is no longer verified. 



We have integrated numerically the system of equations ([16]) for the same param- 
eter values of Fig. [S] and the same initial conditions, but with e = 0.4 and 8 = 0.4. 
The numerical results are presented in Fig. |9] and, we conclude that the anti-phase 
synchronized state still exists for large values of the parameters e and 5. 

Fixing £ to the value e = 0.4, the anti-phase synchronized state still persists for 
5 € [—0.1,3.5]. For 5 = —0.2, and the same initial conditions as in Fig.|9] the system 
in-phase synchronizes. 

The comparison between figures Fig. |9^) and Fig. ^) shows that the definition 
of anti-phase synchronization used in the previous section is specific to the case of 
non-linear oscillators with equal parameters. In Fig.|9]3), the two pendulums clearly 
synchronize in anti-phase, but, lim,„^oo 0i(fn) ^ — lim,„^oo 92{t„). 

As the model presented here shows anti-phase and in-phase synchronization phe- 
nomena for oscillators characterized by different parameters, these numerical results 
show that the equality between the eigen-periods of the two pendulums is not re- 
quired to obtain synchronization. 
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We have proposed a model of interaction between oscillators leading to exact anti- 
phase synchronization. This phenomena has been observed by the first time, in 1665, 
by Christiaan Huygens. 

The interactions parameters in our model are a damping constant p, a stiffness 
constant ^ of a linear string, and a mass parameter M. This mass parameter is as- 
sociated with the interaction, and not with the individual oscillators. For moderated 
values of the (wet) damping constant p, k> 0, and M > 0, the asymptotic solutions 
of the model equations converge to a stable limit cycle in the eight-dimensional 
phase space of the model equations. Numerically, this limit cycle is the only stable 
attractor of the dynamics of the interacting oscillators. 

For smaller values of the damping constant p, two stable limit cycles in phase 
space coexist. One corresponds to the anti-phase synchronized state of the two pen- 
dulum clocks, and the other corresponds to the in-phase synchronized state. The two 
limits cycles are reached by different initial conditions in the model equations. 

The transition between the one-limit cycle solution and the two-limit cycle solu- 
tion appears by a global bifurcation tuned by p. After the bifurcation, the anti -phase 
state coexists with an approximate quasi-periodic in-phase state. 

Changing the parameters of the individual pendulum clocks, we obtain the same 
synchronization properties as in the case of oscillators with identical parameters. 
This fact shows that the interaction mechanism purposed here is robust to changes 
in the parameters of the nonlinear oscillators. 

In all the cases analyzed numerically, the anti-phase and the in-phase synchrony 
occurs with periods different from the eigen-periods of the individual oscillators. 
This shows that the equality between the eigen-periods of the individual oscillators 
is not required to obtain anti-phase or in-phase synchronization. 

An important new issue introduced in the model is the possibility of existence 
of small movements of the attachment points of the pendulums clocks, a situation 
clearly avoid in the modern experimental devices, and diverging from mechanism 
of synchrony proposed by Kortweg, |9|. This explains why modern experimental 
setups have not been able to reproduce the original Huygens's results. 
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